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METHOD AND APPARATUS FOR NOISE FLOOR ESTIMATION 
Field of the Invention 

The invention relates to the field of telecorainiiJiications and in particular the si^al 
processing in receivers of coramunication signals, 

5 Background 

In the art of telecoinmimicatioiis, it is desirable to estima.te the average noise power 
contamed in a received signal. Typically, the received signal is a frequency division 
multiplexed composite signal containing several distinct signals. It is important to be able to 
estimate the noise power without a priori knowledge of the number, nature, carrier-frequency, 
bandwidth or amplitudes of the distinct signals. 

In the art^ noise is considered an additive phenomenon. It distorts a signal by adding 
random time- varying values to the amplitude qf the uQia-uoi3e component of the signal. It is 
thus of great importance to develop a radio receiver that is robust to the noise phenomenon. 

The general problem of estimating the variance of the additive noise in a communication 
15 signal composed of many narrowband frequency division multiplex signals is important in 
many radio coramimication systems. It allows setting the detectors threshold to its maximum 
setisitivity while limiting the probaibility of a false alarm to a value legs than or equal to a 
desired level. The detector decides that a signal is present when the energy of a signal, over a 
period of time, exceeds a threshold value. The detected signals can then be passed to a 
20 receiver or to measurement equipment, for demodulation or to acquirer more information 
about the signal. 

Currently there are several techniques that are used in diis pursuit, but all come with 
sacrifices. One such technique is to reserve a portion of the bandwidth of the composite 
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signal to carry a "noise channel". Thus by measuiing the power in this reserved baad, an. 
estimate of The noise power in the other channels can be obtained. It is coxnmQA though for 
these "Boise channels" to be difficult ro identify as their placement in the spectnim of the 
composite signal is nst governed by any global convention which leads to different channel 
assignments based on geograpliical location of the measurement equipment Other errors are 
introduced when the "noise channel" is affected by the sigB^ls in adjacent channels. Tills 
results in a biased estimation of the noise floor If the channel assignment is unknown, it may 
not be possible to identify the ^*noise channel" and determine the noise floor power leveL In 
addition to these problems the noise floor tends to fluctuate with environmental and operation 
conditions. 

In the artj it is cammon practice to set a noise tbxeshoJd value, determined in large part by 
the estimated noise floor, whereby a signal is considered present if and only if it is above the 
noi^e threshold. 

The traditional approach is to set a threshold value depending on the user needs. YvTien 
high sensitivity is required, the ttireshold is set to a valae close to the user observed noise 
floor level. When only strong signals are of interest, the threshold is set to a high value 
relative to the observed noise floor leveL This approach is highly subjective and may result in 
being too aggressive or too conservative for the actually measured signal environment . The 
threshold is typically set manually in a commmiications system. This is not feasible when 
fast response tq dynamically changing conditions is required. Very few automated 
techniques have been studied to address the general problem of rapidly and dynamically 
estimating the noise floor level or the signal-to-noise ratio. One typical method is to isolate a 
channel and to assign this channel as a noise only channel for estimation purposes. This 
method is difficult to apply in many new digital coimnunication systems where the resources 
are often atilized close to their capacity. It is also a problem for dynamic signal environments 
based on multiple access. Another typical method is tq assume that the composite wideband 
signal non-noise subspace dimension is known. This is again difficult to realize in practical 
situations with unknown signals where robustness is a key feature. Qn another front a 
technique based on mqipha logical binary image processing operators (similar to rank-order 



filtering) as described in "Automatic Noise Ftoor Spectrum Estimation in the Presence of 
Signals", M.L Ready, M,L. Downey, and LJ. Corbalis, Proceedings of the 3r^ Asilomar 
Conference OK Signals, Systems and Comput^rs^ voL 1, pp 877-881, November 1997, 
processes a binary image of the widebajid received power spectrum of a. composite signal has 
5 "been proposed. ThuS;. the approach does dot process the tirae samples or frequency domain 
traiisfqmied data directly, bnt the instead it processes an image of the spectrmrL This involves 
two dimensional signal processing that is usually complex. Furthcrxnore, practical 
performance results h^-ve yet to be quantified with such an approach- In general, techniques 
based on textbook results tend to assume several conditions that are not truly present an many 
10 bands of the radio spectrum. 

The computational complexity of a method is a limiting factor of the method's 
applicability^ Typically, a spectrum jnonitoring receiver cannot observe the full spectrum it 
is allotted, thus it divides the spectrum into smaller blocks that can be scanned, each of these 
blocks possibly contains a plurality of signals. Each block is received sequentiaU3'^, an 

1 5 estimate of the noise JQoor is arrived at, and then the receiver shifts its frequencj'- to receive 
another block. This cycling through the allotted bandwidth must be done rapidly to ensure 
that all messages being transmitted in the different parts of the allotted spectrum can be 
received. If the method to derive a noise floor estimate is too complex then it either cannot 
keep pace with the shifting of the block or it causes the shifting of the blocks to slow down, 

20 in the event of either result the method ceases to be use&l. 

Thus, it would be desirable to design a system that resolves the problems that tlie above 
mentioned methods have with narrowband signals of unknown nature^ as described in ^TSToise 
Floor Estimation for Wide-Band FFT Filter Banks". It would also be desirable for a nev^' 
system to have fliUy digital characteristics so as to permit real-time signal performance^ 
25 without requiring the dedication of a reserved firequency band for use as a noise channei. 
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Summary of the Invention 

It is ^ object of the invention to remedy numerous problems found in the prior art 
5 methods. 

It is anotbear object of the idLveation to provide a noise floor estimation method and 
apparatus capable of operation on a wideband signal without fiiU knowledge of the nature of 
th& narrowband signals contained therein. 

It is another object of the invention to provide digital characteristics to allow real-time 
1 0 signal processing v^fli low levels of computational compledty. 

The preferred embodiment of the method of the invention, with very little a priori 
knowledge of the signals^ allows both robust perfonnaace in non-ideal environments, and an 
automated maximuiu sensitivity threshold setting method. These modifications make the 
process of the invention an attractive solution to problems emerging in the field of digital 
15 conamunications. As an added benefit, information derived by the method of the invention 
may be used in spectrum quality indication, interference analysis^ and signal demodnlation. 

In the preferred embodiment of the method in accordance with the present invention, the 
noise floor is estimated by carr^TOg cur the steps of: 

a. representing the continuous signal as a series of discrete frequency and 
2 0 amplitude values; 

b, creating a histogram based on the discrete frequency and amplitude values; 
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c- dfiri\ang a noise floor estimate from characteristics of the histogram. 

The discrete frequency and amplitude values are preferably obtained by the step$ of: 

d~ sampling a received wideband sign?il by a plurality of analogue-to-digital 
convertears to generate a series of output signals; 

5 e. windo^'ing the output signals of the analo^ue-to-digital converters; 

f. applying a mathematical tr:3nsfomij preferably a Fast Fourier Traosformj to the 
results of the windowing to ohtain a series of discxote frequency values; 

coarerting an amplitude of each discrete frequency value to log-domain 

representatiorL; and 

"I 

10 h- rounding the log-domain representation of the amplitude for each di$crete 

frequency value to the nearest integer valu6 to generate a discarete amplitude 
value. 

The "Windowing step preferably includes the steps of. 

i. selecting a discrete weighting function; 

15 j , multiplying the value of each output signal of the series fay a corresponding 

element of the discrete weighting function. 

The amplitude of each value of the discrete frequency series is preferably coTiverted to log 
domain representation by multiplying 20 by the base 10 logarithm of the xnagnitude of the 
element; whereby the log domain representation of the amplitudes results in the amplitudes 
20 being escpressed either as decibel (dB) values; ox as decibel milfrwatt (dBm) values. 
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In a first preferred embodiment of the steps of creating and analysing the In'stogram to 
estimate the noise floor, herein referred to as the discrete level crossing rate method, a 
histogram chart is constnicted where the minimuiti v^lue is equal to the minim-urn integer 
value in dB of the received signal strength, and the maxirniun value is equal to one less than 
5 the highest integer value in dB of the received signal strength. Every time a segment of the 
signal is found to have a positive slope, the histogram bins crossed by the segment are 
incremented by one. TTins, for example, if a given segment starts at -1 10 dB and rises to -105 
dB in ttte next sample, die values of the bins representing -1 10, -109, -108,. -107 and -106 dB 
are incremented by one. After the complete histogram has b^en assembled, the next step is 

10 preferably to determine the location on the dB scale of the noise floor level. Typically, the 
histogram will be bimodal in the presence of noise and a plurality of signals, with tbe first 
true local maximum indicating the noise floor. To find the noise floor^ preferably the first 
local maximum starting fi-om the minimum quantised dB value is selected. Because the 
histogram is created from sampled data affected by a random source it is possible to have . 

15 nnmerous local maxima that do not reflect the trae noise floor, as a resnlt it i& preferred to 
check if the local maximum is in &ct the noise floor- In the most preferred embodiment, the 
metliod for the maximum value fonnd includes the step of checldng if a histogram value 
witiiin y dB is higher than the current maximum, where Y is a prefeired number derived 
empirically fi-om the characteristics of the system, most preferably between 2 and 3 dB. If 

20 the no value in the ruext Y db is greater than the value at the maximum then the current 
maximum abscissa is indicative of the noise floor level estimate per bins, iiowever if the 
\ralue is greater Y db away firom the maximum than it is at the maximum then it is cleat that 
the noise floor has not been found. In the event that the maximum is not the noise floor^ the 
next local maximum is preferably used and the Y dB test applied again until the test fails to 

25 produce a value higher than the value at tlie maxifflum. The obtained abscissa value is finally 
preferably corrected by a fixed offset of X dB, where X is preferably dependent on the 
selected window, to give the noise floor level estimate per bin for the processed block of data. 



In a second preferred embodiment of the steps of creating and analyzing the histogram to 
estimate the noise floor, herein referred to as the penalized quasi log-likelihood method^ a 
30 histogram of quantised frequency domain amplitudes is constructed with each bin 1 dB in 
size- The histogram is then transformed into a sorted vector of linear amplitude elements in 
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decreasing order. The sorted linear vector is then preferably reduced firom the size MK to M 
by summing a grGup of K consecutive samples of the sorted vector, la. this method M is 
typically tnxich. larger than K, A quasi log-likelihood fimction L(k) is thea preferably applied 
on the sorted group of M elaments. The negative of the quasi log-likelihood function i$ 
5 preferably added to a pcnalt>' polynomial ftmction p(k) giving an expressiou of the form 
- L{k) ^Cp[k), where C is a constant, preferably empirically determined and most 
preferably near -2.6. The mimmam of ftie -L(k)^ C p(k) function is, preferabljr selected 
where the index k of the minimum amplitude is denoted by q^. If a line is dra^Ti on the 
histogram at the index , it can be seen that the received signal can he di\rided into 2 
1 0 groups^ the first being the message subspace, where the information contained in the non- 
noise sign^i exists J and the noise subspace, where the random noise values exist- Finally, the 
noise floor level estimate per frequency bin is prefmbly obtained by dividing by K the 
average M -<ff,jr -1 smallest values of theM sorted group. Thus the noise floor estimate per 
K firequDcy vin is the average value of the power along the dimensions of noise subspace, 

15 The application where the proposed methods have been used is in spectrum monitoring of 
communication radio signals. The methods arc particularly suitable where fast ^canning 
operations are necessary and an adaptive threshold value is rec^uired. They can also be used to 
automate lengthy manual noise level measm-ement procedures, and the resulting noise floor 
level estimates can be used as an indication of the spectrum quahty. 

20 As stated earlier, it is an object of the invention to derive an estimate in real-time of the 
noise floor level in a wideband signal preferably without laiowiedge of the signal content in 
terms of carrier frequency, bandwidth, md power, and using only a short^time signal. Lathe 
preferred apparatus ia accordance with ll\e inventions the wideband signal, containing a 
plurality of narrowband signals, is preferably received from the analogue output of a radio 

25 receiver and sent to a digitizer where a block of output data is gmerated over a lame period. 
The sampled amplitudeSj thoagh discrete, are real valued and thus the signal is still 
considered continuous in amplitude, Tlie block of output data is then windowed in a 
windowing module and sent to a ohannelizer where the wndeband signal is representied by a 
collection of complex amplitudes at several discrete frequencies. The amplitudes of the 
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signals at tlie discrete frequencies are converted in a converter module to the log dom?an 
representation (decibels (dB) or decibel milliwatts (dBm)) and rounded to ih^ nearest integer 
values. Thus the received signal Is now represented by boih discrete amplitude and 
frequency values. 

5 Brief DoscriptSon of the Drawings 

Exemplary embodiments of the invention will now be furtheoc described with 
references to tlie drawings in which same reference nmnstials designate similar parts 
throughout the figures tliereof^ and wherein: 

Figure 1 is a block diagram of the signal pmccssing elements for the discrete level 
1 0 crossing rate (DLCR) embodiment of the method according to the invention; 

Figure 2 is an example of the discrete level crossing rate with a noise otdy input 

signal; 

Figure 3 is an example of the discrete level crossing rate with both noise and non- 
noise sisals at the xeceiver's input 

15 Figure 4 is an example of the behaviour of tlie DLCR method as a fimction of the 

occupancy; 

Figure 5 is an example of experimental results showing the time bebavionr with low 
and high signal frequency band occupancy for the DLCR method in the 800 MHz band; 

Figure 6 is a block diagram of the signal processing elements for the penalized qua^i 
20 log-likelihood (PLLM) embodiment of the method according to the invention; 
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Figure 7 is a plot of the quasi log-lilceUhood function and the associated penaliziiis 
polynomial; 

Figure 8 is an example of the PLLM(k) fractian for a nosie only aigxial, and a non- 
noise with noise signal; 

5 Figure 9 is an example of the belia-viour of the PLLM method as a function of the 

f-equency band occupancy; 

Figure 10 is an example gf experimeDtal rfesults showing the time beha^dour with low 
and high signal occupancy for the PLLM method in th« SOO MHz band; and 

Table 1 shows the offset for the DLCR method 

10 Detailed Description of the Preferred Embodiment 

In the context of spectrum monitoring^ the input signal contains several channels ttet 
are observed fay a wide-band chaimelisation receiver. The receiver requires kno^'ledge of the 
noise fLoor^ m order to set a threshold that will produce a constant detection false alarm rate, 
when the only assumption is an input signal with additive noise nonnally distributed. For 
15 simplicity" and speed, the receives: must estimate the noise flogrr level fron? the current 

sampled uiformation. vector, which is provided by m an^ogue-to-digital (A/D) converter. 

The time over which the input vector is sampled can vaiy depending upon the 
resolution required, the signal processing stages that follow and the physical constraints of 
moving the data to the subsequent stage. To ensure a fast scanning opearation of the 
20 channelised receiver, it is impoxtant not to use more samples to extract the necessary 
infomiation flxsm the input vector for the noise floor level estunation process than the 
channelised receiver needs. The presence and parameters of these non-noise signals may not 
be known by the monitoring equipment. The a priori knowledge of the receiver is jninimal 
and is lirnited to the information that signals are not wider than the sampling band^^dth. 
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Fqx a first embodiment of the present method, herein refsrred to as the Discrete Level 
Crossing Rat^ method, the signal path of the estimation methgd is presented in Figure 1 . The 
signal to "be analyzed 101 has prieviously been JSltered by methods known in the art to prevent 
aliasing of the fiequency batids after the sampling process. The input analogue signal 101 
5 contains additive noise emmiating from various sources and non-ncyise signals. A plurality of 
amlqgue-to-digiial converters 102 digitizes the signal and stores a series of samples in a 
buffer represented by the vector r 105. The vector r 103 is then windowed or polyphase 
filtered to produce another vector to he sent to a digital filter bank 104. The window choice 
may depend on specific requirements relative to a particular signal fenvironment or on filter 

10 bank design criteria. It is also important th^t any gain m the signal path (noise and non-noise 
sign^s) up to the filter bank output be identifie4 cancelled or compensated, so that the output 
of the filter bank is representing the tme signal power at the antenna input. It is therefore 
necessary to ensure that the equipment is properly calibrated to avoid any biased estimation 
of the noise floor leveL Jn the digital domain, the filter bank is often implemented as a Fast 

15 Fourier Tiansfonn (FFT) to minimise computational complexity, but otlier transforms are 
also possible. The FFT is the prefenred filter bank in this inveution. The filter bank output 
represented by the vector z 105 is then transformed in a transforming module 106 to the 
integer log magnitude domain (dB or dBm) with a rounding operation on 20 times the 
logarithm base 10 of the magnitude of each complex: point of the vector z. The 

20 implementation of the transformation step 106 is simplified by the rounding operation such 
that the precision required in the log computafon can be small. The resulthng vector x 107 is 
then processed by the Discrete Level Crossing Rate (DLCR) module 108. The analogue 
Level Crossing Rate (LCR) fimction is well knov^ hx the signal pmcessing community 
where it is mainly used in wave propagatiqn and cotnmunicatiou channel modelling to obtain 

25 statistical information on the rate of change of feding processes on tran^mitbed si^als^ and 
can be found in munerous standard texts including "Mobile Communications Engineeiing" 
W,C-Y. Lee, McGraw-Hill, 19S2, The LCR h a second-order statistical operation that is 
time-dependent. It is usually sufficient to consider onl}'' the upward or the downward crossing 
in the ervaluatiou of the function. The LCR is defined as the total number of expected 

30 crossings of the signal envelope, at a givai amplitude A, per second. The standard LCR 

operates on time domain signals and is a cDntinuous level function with a contumous signal at 
the input The current signal has a dascrete-fiequency discrete-amplitude representation. To 
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overcome this difficulty, the solution is to use a discrete aniplitude version of the LCR in the 
form of a histogram containing the number of times a discrete decibel amplitude valu^ is 
crossed in the upward direction. Therefore, with discrete-frequency signals, an approximation 
to the true LCR is built as a histogram of the number of oocurrences of the positive crossing 
5 (positive slope) of an integer arupjitude value In dB. For an additive Gaussian noise only 
signal, tixis procedure produces a toiction having a shape of the form , 



n[Z^A] = 



where A is the integer level in dB at wliich the number of crossing is required^ pis a 
correlation parameter of the frequency levels, aad is the variance of the additive complex 
1 0 Gaussian noise process in an PFT bin. The fimction above has a single rnaximum that can 
easily be verified to be located at lOlogj^ (cr" ) . The parameter to estimate is the value of cr^ 
and it cm be esthnated with the position of the maximum of the function nlA] when only 
noise is present 



The situation where one ot several non-noise signals are present with the additive 
1 5 noise needs additional signal processing. For tids case, it is possible that the marimum of the 
ftttiction is not the noise floor level since ftie signal might cross a given amplitude more often 
if the observation bandwidth h highly occupied wth non-noise signals. To avoid potentially 
large bias errors, the first local maximum, starting from the lowest dB values, is selected- To 
check if this local ma?timum is not due to a random fluctuation in the construction of a given 
20 discrete LCR histogram, it is then verified that there is no greater histogram value on the right 
side of the local maximum closer to FdB relative to the current local maximmn- This 
verification also stabilises the method by preventing large overestimates of the noise floor 
level when the channel occupancy is high. The search for the appropriate maximum 
perfcrmed in search module 110, on the output of the Discrete Level Crossing Rate module 
25 108, which is referred to as 109, can have different forms. The value of ^depends on the 
d>'namic range of the non-noise signals. For this application, a value of 3 dB is used. Other 
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tnechanistns to stabilize the method in presence of high gccupancy narrowband sign^s are 
possible. The abscise value of the mg^cimum; found in search module 110 is the noise floor 
level estimate 1 J 1 per FFT hiri within an offset. The offset to be added to the abscise value 
depends on the window or the polyphase filter used. Table 1 gives the offset for four popular 
5 windows. Thus the noise floor level estimate per FFT bin in dB is 

[dB] = position of appiopriate maximmn [dB] + Offset [dB]. 

The noisfe floor level estimates S-~ can be average for the sartie frequency band to 
reduce the variance if the noise floor level is known to vary vei^*" slowly. 

To evaluate the performance of this method, simulation software was employed and a 
10 real-time implementation was developed to gather statistics about the noi$e floor ievel 
estimates. The performance in terms of the bi^ and the variance of the noise floor level 
estimate of the inve^ited method was presented as a function of the channel occupancy^ and as 
a fimction of the information signal power for a scenario of a specific noi?-noise signal. A 
sampled FM signal in the 1 39 to 144 MHz band wa^ used as the transmitted signal in the 
15 simulation software. The channel bandwidth was set to 30 kHz and che windowed FFT filter 
ba-ok used a Blackman window. The transmitted signals all have the same average power at 
the receiver input. 

Figure 2 demonstrates the constniction of the DLCR histogram for a particular vector 
s when only a noise signal is present hi the simulation data at the input. In Figure 2 a, a 

20 Blackman window is used with 1024 filter bank outputs foimded to their nearest integer 
values. The resulting DLCR is preheated in Figure 2b where the maximum 201 of the 
histogram is at dB. Adding the 3 dB offset gives a noise floor level per FFT bin of -6 dB^ 
the true value. Figure 3 demonstrates the construction of the DLCR histogram for a particular 
vector X when a composite signal having a channel occupancy of 12.5% with K== 4hins per 

25 channel is present with the noise signal in the simulation data at the input of the A/D. The 
signals have a signal power 20 dB above the noise floor level. In Figure 3a, aBLacfcman 
window is used with 1024 filter bairic outputs rounded to their nearest integer values. The 
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resulting DLCR is presented in Fig. 3b where the maxiroam of the histogram is now at -S dB 
301. With the 3 dB offset, the noise floor level per FFT bin is of -5 dB, A bias of 1 dB is 
introduced in Ihis example^ In genial, the algorithm will introduce a small bias as a function 
of the narrowhmd non-noise signal spectrmn shape,, the signal amplitude and more 
5 importantly the channel occupmcy. hx ihe following, the number of channels M has been 
fixed to 64 and the nwnber of bins pea: channel to 8, These nmnbers are typical of what is 
expected in practice. Other choices of number of channels and nnmber of bins per channel ^re 
possible with different perfoxmance results- Figure 4 illustrates the perfoimances for the 
^'alnes of signaJ-to-ngise ratios of 10, 20, and 30 dB in 30 kHz with the true value of the noise 

10 floor level being at 0 dE. Each column bar is the average of 1000 independent estimates with 
an error bar for the standard deviation. Figure 4 shows lhat the estimate bias is increasing 
with the channel occupancy. It is within 3 dB of the true value up to 68.75 % of channel 
occupancy 401. Figure 5 presents the real-time impleraesitation results where 1000 estimates 
were computed as a function of time for real-life signals. The signal of interest ws over a 4.5 

1 5 MHz bandwidth sampled at 5.12 MHz with 5 kHz bin resolution. The size of the vector r is 
thus 1024. Two curves are presented: one where the centre fireqnency was at S4S MHz with a 
low occupancy of approximately 3 %, 3nd one where the centre frequency is 880.05 MHz 
where tlie occupancy is about 54 % with an average signal power of 30 dB above the noise 
floor level. Since the two baiids are very close to each other, it is assumed that the noise floor 

20 levels of both are the sanie. The algorithm shows a bias of about 2.2 dB for an occupancy 
level of 54 % relative to the low channel occupancy band. The variance pIso increases when 
the occupancy increases- In this particular case, die variance is tripled. The performance of 
the method of the invention is thus very good considering that no a priori knowledge of the 
signals is available to the receiver. 

25 Another embodiment of the noise floor level estimation method of the invention is the 

method presented in Figure 6, which is referred to herein as the penalized log-likelihood 
method- The first seven steps of the method (601, 602, 603, 604, 605> 606^ 607) are the same 
as the first seven steps (101, 102, 103, 104, 105. 106, 107) of the DLCR method of the 
embodiment presented in Figure 1, Thns, the received sigtial 601 to be analyzed is assumed 

30 to have been filtered to prevent aliasing of the frequency bands after the sampling process. A 
phirality of analogue-to-digital converters 602 digitizes the signal and stores a series of 
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samples 'm a buffer represented by tlie vector r 603, The vector r is then windowed or 
polyphase filtered before being applied to a filter hmk 604, It ig also important that any gain 
in the signal path (noise and non-noisc signals) up to the filter hank output be cancelled or 
compensated^ so that the output of the filtet bank is representing the true signal power at the 
antenna mput It is therefore neoessary to ensore that the equipment is properly calibrated to 
avoid any biased estimation of the noise floor level. In the digital domain, the filter bank is 
often implemented as a Fast Fourier Transform (FFT) because of its cowiputational 
efficiency. The filter bank output represented by the vector z 605 is then converted to the 
integer log magnitnde domain 606 with a rounding operation on the product of 20 times the 
logarithm base 10 of the magnitude of eaeh coxnplex point of the vector z 605, Here again, a 
requirement is that the noise power level be detenniQed from a single FFT trace. The rounded 
dB values in the vector x 607 are then used to build a standard histogram 608 with histogram 
bins corresponding to the integer dB values. The histogram computed is fhe^i used to build an 
ordered sequence 608 of amplitude dB values in decreasing order w 609, The vector is then 
of the same len^ and contains the same elements as the vector x but with elements in 
decreasing order. Thus, the block of histogram and the decreasing ordered sequence 608 is 
equivalent to a decreasing sort operation on the rounded dB values. The present approach is 
much less complex in terais of number of operations required- to perform the sorting than a 
true sorting method^ The decreasing sorted integer dB values are then transformed back to the 
linear dotnaiB 610 with the transformation 



for each of the elements in the vector w 609. The decreasing sorted linear values are then 
summed by groups of 610 starting ^vith the largest values. This grouping reduces the 
dimension of the problem from iWK" to M. Another impact of the grouping operation is the 
tendency of the group sample power to he de-cptrelated even mth the use of a window 604^ 
because the decreasing ordered samples are not necessary from adjacent FFT bins. The 
following block 612 is th^ computation of the penalized quasi-log-lifcelihood function with 
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s 



the expression 



where is a quasi Igg-likelihood function of the M sorted group elements, C is a constant 
and p{K) is apeixalty fiuiction. In tins embodim^t, the quasi-log-likelihood function is given 
5 by 



where U are the ^decreasing valaes of the power in the grouped iTbins 611^ thou^ it should 
be noted that a variety of log-likelihood, and quasi-log-likelihood functions could provide 



10 variance of independent and zero me^n Gaassian random variables^ but the elements of a 

received signal are correJated and thus are not independent nor are fliey necessarily zero mean 
Gaussian random variables. The penalty function /?(ft) is added to the negative of the quasi 
log-likelSiood fimctioii to form the function PLLM(A) over which the minimum T.vill be 
searched 614. The penalty function p(fr) is often related to an Information Theoretic Criteria 

15 as described in *tD etection of Signals by information Theoretic Critena'' M- Wax and T- 

KfliJath, IEEE Transactions Aconstic, Speechj Signal Processing volume 33, Dumber 2, pages 
387-392, April 1985, The penalty ftmction in this embodiment, however, \% different £bom 
those derived earlier in tltat the coef jScients U are not truly statistically independent, and that 
the intent is not to determine the foU model order as in the reference cited above. For this 

20 embodiment, the penalty ftmction is a second order polynomial given by 




acceptable results. The function would be a true log-likelihood function if //-were the 
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other polynoniials or functioiis are possible, but this second order polynomial is adequate for 
a good performance and complexity trade off. The polynomial function is scaled by a 
constant C to avoid erroneous detection due to the variance in the data. The value of C has 
5 been detenxiined empirically and may vary depending on the final goal of the design or the 
signal environment at the input. For this application, C — -2.6 is the preferred value. The only 
assumption to determine tlie above polynomial is that the FFT coefiScients are distributed 
according to a Gaussian probability density function. If the additive noise time data samples 
do not have a Gaussian distribution the linear combination of the FFT will tend to trmjsforai 
1 0 them into variables having a probability distribution close to Gaussian v^ables for practical 
and large dze of MK- The method is therefore robust to non-ideal conditions where the 
additive noise might not be truly Gaussian, The scale factor C gives a degree of freedom in 
the optimization process. The output of function PLLM(A;) is then searched in step 614 for the 
index of its mnimum defined as 

15 ==argmin{-i(^)4-Q;fe)}. 

The above expression means that q^^ is the index at the mimmum of the function 
?l-LM{k), The noise floor astmiate per A' FFT bins; is the average of fhe M-q.f^ ~1 smallest 
numbers of the sorted vector values /,-. The noiso floor estimate 6 15 per FFT bin is then given 
by 



J 



I M-1 



Tlie noise floor level estimates can be averaged for the same frequency band to 
reduce the variance if the noise floor level is known to vary very slowly. 
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To evsQuat^! the performance of the PLLM algorithm, simulation software was employed 
and a real-time LmplenQentation was developed to gather statistics about the noise floor level 
estimates. The performance in tenns of the bias and the variaiice of ihe laoise floor level 
estimate of the taventfid method is presented as a flmction of the channel occupmcy and as a 
5 function (if the infoimatiou signal power for a scenario of a specific signal. A ssanpled FM 
signal in the 139 to 144 MHz band is used as the transmitted signal ici the simulation 
software. The channel bandwidth is set to 30 kHz and the windowed FFT filter bank used a 
Blackman window. The transmitted signals all have the same average power at the receiver 
input 

\0 Figure 7 illustrates the tiormalized quasi log-likeUhood fimctioa 701 aad the second-order 
polynomial approximation 702 with the scale factor C - 1 for ^ Gaussian noise only input 
signal. As the figure shows, the two curves are similar over the range of the index k Figure S 
presents the PLLM(^) function for a noise only, and a non-noise plus noise signal. It can be 
seen that the index of the minimum of the PLLMCJt) function is indicative of the signal 

15 subspace dimension. Since the remaining subspace includes noise only components, then Ifae 
average is an estimate of the noise floor level. In the following, to reduce tlie amount of 
simulation results to present, the number of chamiels If has been fixed to 64 and the number 
of bins per channel K to 8. These numbers are typical of what is expected in practice. Other 
choices pf number of channels and number of bins per channel are possible with different 

20 performance results- Figure 9 Illustrates the performances for the values of signal-tc-noise 
ratios of 10, 20, and 30 dB in 30 kHz with the true value of the noise floor level being at 
0 dB. Each columo bar is the average of 1 000 independent estimates witii an error bar 
representing the standard deviation- Figure 9 shows that the estimate bias increases with the 
charmel occupancy. The estimate is within 2 dB of the tme value up to 84.375 % of channel 

25 occupancy 901, Figure 10 presents real-time implementatioa results where lOOO estimates of 
the noise floor level were computed using real-life signals. The signal of interest was over a 
4,5 MHz bandAvidth sampled at 5.12 MHz with 5 Miz bin resolution. Two curves are 
presented: one where tlie centre frequency was at 848 MHz with alow occupancy of 
approximately 3 and one where the centre J&equency is 880.05 MHz where the occupancy 

30 is about 54 % with an average signal power of 30 dB above the noise floor leveL Since the 
two bands are very close to each other, it is assumed that the noise floor Ie\'els of both are the 



17 



same- The algorithin shows a bias of about 1 . 15 dB for sn occupancy level of 54 % relative to 
the low channel occupancy band. The variance also increases when the occupancy increases- 
In this particular case, the variance is almost tripled. The performance method is thus very 
good considering that no a priori knowledge of the signals is available. 



5 The exemplary embodiments described above are provided for the sake of clarification 

and are not intended to limit the scope of the invention. Of coxirse^ numerous variations and 
adaptations may be made to the abovs-described embodiments of the invention, without 
departing &om the spirit of the invention^ which is defined in the claims. 

Chmiges and modifieations in the specifically deseribed embodiments can be carried out 
1 0 without departing from the $cope of the invention which is intended to be limited only by the 
scope of the appended claims. 
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WHAT IS CLAIMED IS : 

1 . A method of estimatuig the noise floor of a wideband analogue signal comprising 
the steps: 

a. representing the continuous signal as a series of discrete frequency and 
amplitude values; 

creatirxg a histogram based on the discrete frequency and amplitude values; 
c. deriving a noise floor estimate from characteristics of the histogram. 

2. The noise floor estimation method of claim 1 wherem th^ meihod of obtaining the 
series of discrete frequency and amplitude values includes the steps of: 

a. sampling a received wideband signal by a plnrality of analogue-to-digital 
converters to generate a series of output signals; 

b. Nvindowitig the output signals of the analogue-to-digit^l converters; 

c. applying a mathematical tratisforra to the results of the windowing to obtain a 
series of discrete frequency values; 

d. converting an ampBtudfi of each discrete frequency value to log-domain 
representation; and 
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5. roundrng the log-domain representation of the amplitude for each discrete 
frequency value to the nearest integer value to generate a discrete amplitude 
value. 

The noise floor estimation method of claim 2 whereiu the windowing process 
includes the steps of: 

a. selecting s discrete valued weighting function; 

b, multiplying the Y^ue of each output signal of the series by a cofresponding 
eleiDfint of the discrete weighting function. 

The noise floor estimation metiiod of claim 2 whearein the jnattiematic&l traxisfonn 
used is a Fast Fourier Tiaiisform. 

The noise tloor estiiuatioti method of clatni 2 wherein the amplitude of each value 
of the discrete frecjuency series is converted to log domain representatioTi by 
multiplying 20 by the base 10 logarithm of the magnitude of the element 

The noise floor estimation method of claim 2 wherein the log domain 
representation of the amplitudes results in the amplitudes being expressed as 
decibel (dB) values. 

The noise floor estimation method of claim 2 wherein the log domain 
representation of the amplimdes results in the amplitudes being expre^ed as 
decibel milliwatt (dBm) values. 
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The noise floor estimation method of claim 1 whereiu tlie method of buiWing the 
histogram based on the discrete frequency and amplitude values includes the steps 
of: 



a. establishing a lowest bin representixig the lowest integer dE value of the 
discrete series representing the wideband signal; 

b. ostahlishiBg a highest biti repiesenting the highest integer dB value of the 
discrete series representing tiie wideband signal; 

c. establishing bins for each integer dB value between the lowest and highest 
bins; 

d. incrementing the value of each bin when a segment of the series f epresentibag 
the wideband sigD.al crosses the bin with a positive slope. 

The noise floor estimation method of claim I wherein the stop of deriving the 
noise floor estimate from the characteristics of the histogram includes tlie stops of: 

a- defining the lowest dB bin as a starting point 

b. determining the next lowest valued local maximum on the histogram; 

c. performmg a Y test; 

d- repeating steps h aiui c until the Y test fails^ 

e. setting the noise floor by adding an offeet to the dB value of the msxiraima of 
the histogram that caused the Y test failure. 
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10. 



The noise floor estimatioTi method of claim 9 wherein the Y test includes the steps 
of: 



a. exaraining all poiats in the next Y dB ; 

b. considering the test a pass when a point exists in the next Y dB which has a 
5 higher v^lue than the starting point; 

c- considering the test a fail when no point exists in the next Y dB which has a 
higher vain? than (ike starting point. 

11, The noise floor estimation method gf claim 10, wherein Y is 3 dB. 

12, The noise floor estimation method of claim 9, wherein the offeet is detennined 
10 based on observed characteristics of the signal and the windowmg pr(?cess' 

discrete weighting function. 

13, The noise floor estimation method of claim 12, wherein ttie offset is selected 
&oin the group of 2 dB for a rectan^ar window, 2.75 dB for a Hannigg window, 
3 dB for a Blaekman window, and 3 J. dB for a flat top window, 

15 14. The noise floor estimation method of claim 1 , wherein the step of creating the 

histogram includes the steps of: 

a. establishing a lowest bin representing the lowest integer dB value of the 
discrete series representing the wdehand signal; 

b. establishing a highest fain representing the highest integer dB value of the 
20 discrete series representing the wicLeband signal; 
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c. establishing bins for each integer dB value between the lowest and highest 
bins so tet there are a total of MK bins; 

d. increnientiiig the bins for each time an element in the discrete series fells into 
the bin. 

The noise floor estimation inethod of claim 14, wherein the step of deriving the 
noise floor estimate from the characteristics of the histogram includes the steps 
of: 

a. sorting the elements of the histogram in decreasing order of amplitude to 
create a sorted vector; 

b. reduemg the size of the sorted linear vector from MK to M by smrmmig 
^oups of K consecutive elements of the sorted linear vector for achieving a 
more discretised amplitude representation; 

c. applying one of a log-likelihood fimctioiL, and a qmsi log-likelihood jfunction^ 
to the M clemeaats of the sorted hnear vector to achieve a discrete functioti 

d. subtracting L(k) from a multiple (£7) of a discrete penalty function to 
obtain the fimction -LlkJ-i-C p{k); 

identifymg the indesx, denoted by , at which the minimum of the 
- L{k) + C p{k) equation is achieved; and 

L computing the noi^e floor level estiuiate per FFT bin by dividing the mem of 
the M — q^Tf, ~ 1 smallest values of the M sorted vector by K. 



16. 



Hie noise floor estimatioa method of claim 15, wherein M is considerably larger 
than K, 



17. The noise floor estimatioii meHiod of claim 16, wherein M = 64, 



1 S , The noise floor estimation method of clahn 1 6 ^ wherein K=8 , 
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The noise floor estimauon method of clahn 15 wherein is represented by the 



quasi-log-likelihood function L{k)= Khi 



1 



, where k 15 the 



index of the fimction. 



20, The noise floor estimation method of claim 15> wherein the penalty function is a 
polynomial. 



10 2T The noise floor estimation method of claim 1 5, wherem the penalty function is 

represented by the second order polynomial jflmcdon 



{ M-l J [, M-l ) 



3,76 



22, The noise JEloor estimation method of claim 15, wherehi lie constant C is -2,6. 



23.. A wideband analogue signal noise floor estimation apparatus comprising: 



15 



a. a digitizer module for creating a representation of the continwns signal 
comprised of discrete frequency and amplitude valnes; 
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b. a histogram module for generating a histogram based on the discreie fi:equeucy 
and amplitude valines; 

c. an estimatioxi module for deriving an estimate of the Mise floor of the 
wideband signal based on the characteristics of the histogram. 

The noise floor estimation apparatus of claim 23, whea-ein the digitizer module 
fuither comprises: 

a. a stoiplitig module including a plurality of analogue-to-digital converters for 
generating a series? of output signals; 

a windowing module for weighing the output signals of the sampling element 
to generate weighed output signals; 

c, a transfotming module for applying a mathematical transform to the weighed 
output signals to create a si^al comprised of discrete frequency values that 
represent the original signal; 

d, an amplimde domain converter for converting the linear amplitude values to 
Log-domain representation; and 

e, an amplitude disGtetizing module for representing the output of the amplitude 
damain conversion element ^ a sequence of integer valued amplitude levels. 

The noise floor esthnation apparatus of claim 24, wherein the windowing module 
further includes a weightiag element for multiplying oach value of the output 
series by a corresponding element of a preselected discrete valued weighting 
function. 



The noise floor estimation apparatus of claim 24, wherein the transforming 
module applies a Fast Fourier Transform, 

The noise floor estimation apparatus of claim 24, whereii) fhe amplitude 
discretizing module is constructed to convert each amplilude value of the discrete 
frequency series to 20 times base 10 logarithm of the magnitude of the value. 

The noise floor estimation apparatus of claim 24, wherein the amplitude domain 
converts outputs amphtude values as decihel (dB) values. 

The noise jfloor estimation apparatus of claim 24, wherein the amplitude domain 
converter outputs amplitude values as decibel milMwatt (dBm) values. 

The noise floor estimation apparatus of claim 23, wherein the histogram module 
further includes 

a. a low bin establishing element for oreatin^ a low valued bin to represent the 
lowest inieger dB valvj? of the discrete series representing the wideband 
signal- 

b. a high bin establishing clement for creating a high valued bin to represent the 
highest integer dB value of the discrete series representing the wideband 
signal; 

c. a tertiary bin creation element for creating a bin for each integer dB v^ue 
between the lowest and highest bins; and 

d. a bin count incrementing element for incrementing the value of each bin when 
a segment of the series representing the wideband signal crosses the bin with a 
posihve slope. 



The noise floor estimation apparatus of claim 23, wherein the estimation morale 
further includes: 

a, a maxima fiading element fox finding the next left most maximum from a 
given starting point, that in the absence of previous data takes the lowest dB 
bin as a starting pomt- 

b - a Y tesi element for performing a Y test; 

c. a decision element for calling upon the maxima finding element until the Y 
test element reports a fail; and 

d, a noise floor settmg element for providing a noise floor estimate by adding an 
offeet to the dB value repprted by the maxima finding element that caused the 
Y test element to report a fail. 

The noise floor estimation apparatus of claim 3 1 , -whearein the Y test element 
farther mcludes: 

a. an examination element for searching th6 Y dB to the right of the given 
starting point for a value higher than the starting point; and 

b. a reporting element for reporting a fail when vo point exists in the next Y dB 
that has a higher value than the starting point and report? a pass if fhere is a 
value in the next Y dB that is greater in value than the starting peine. 

The noise floor esti^nation apparatus of claim 32, wherein Y is set at 3 dB. 
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34. 



The noise floor estitnatioii apparatus of claim 31, wlierein the offset used by the 
noise floor setting element is based an observed characteristics of the signal and 
the windowing process" discrete weighting functioiiv 



35. The UQi^e Aqox estimation apparatus of Dlaim 34, wherein the offeet is selected 
5 from the group of 2 dB for a rectangular window, 2.75 dB for a Harming window^ 

3 dB for a Blackmaa window, and 3.2 dB for a flat top window. 

36- The noise floor estimation jipparatus of claim 23, wherein the histogram module 
fiirther includes: 

a. a low hin establishing element for creatirvg a low valued bin to represent the 
10 lowest integer dB value of the discrete series representing the wideband 

signal; 

b. a high bin establishing element for creating a fiigh valued bin to represent the 
highest integer dE value of the discrete series representing the wideband 
signal; 

15 c. a tertiary bin creation element for creating bins for each integer dE value 

between the lowest and highest bins; and 

d. a bin count increanenting elem^t for incrementi-ng the value of a bin for each 
time an element in the discrete series falls into the bin. 

37. The noi5e floor estimation apparatus of claim 23^ wherein the estimation module 
20 furtiT.er includes: 



a. a sorting element for creating a vector containing the discrete amplitudes of 
the input signal in decreasing order; 
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b a vector si2:e feducing element for reducing the size of the sorted linear vector 
from MK elements to M elements by smnming grotjps of K consecutive 
elements of the sorted linear vector tg achieve a more discretised amplitude 
representatioD; 

c. a log-likeliliood element for applying a log-likclihood, or a quasi log- 

likelihood fianction, to the M elemeats of the sorted linear vector output from 
the vector reducing element to achiefve a discrete fimction 

d- a penalty function clement for subtracting the discrete function L(k) firom a 
multiple (C) of a discrete penalty function p(k) to obtain the function 
- Z(i)+ C p{k) (PLLM froctlon); 

e. an index id^ntilication element for identifying the indsx at which the 

minimum of the PLLM function, - L{k)+ C p{k)^ is achieved and identifying 
the indeix, denoted by q^^ , at which the minimum of the - L(k)-\- C pik) 
equation is achieved; and 

£ a noise floor setting element for providing a noise floor estimate by di%riding 
the mean of the M — q^^—l smallest values of tlie M sorted ^'octor by K. 

The noise floor estimation apparatus of claim 37, wherein the penalty fonction 
element is constructed to apply a polynomial as the penalty fonctiotL 

The noi$e floor estimation apparatus of claim 37, wherein the peuahy function 
element is constracted to apply the second order polynomial function 




MK as the penalty fiinction. 



ABSTRACT 



A method and device are provided for an improved system and method of estimating the 
noise jaoor of a wideband si^al without full knowledge of the nature of the narrowband 
signals contained thereiB^ and provides digital characteristics io allow real-time signal 
processing with low levels of compntatiotial complexity. The method and device require 
very little a priori fcoowledge of the nature of the sigiials, allow both robust performance in 
uon-ideal enviroixcnentSj and automated maximom sensiti'^'itj^ threshold setting 



30 



0 



3.0xia 



10 



2-5x1(1' 



200 



400 



600 



800 



• Signal and noise (880.05 MHz) 



n£an= l.lE-10 

standard deviati<iii = 3.5E-1 1 



1000 




Noise ODly (84S MHz) 

meaa=6.6E-ll 

standard deviation = 1.3E-11 



I.SxIff 



i.oxia 




400 600 

Time index 



1000 



Fig. 5. 




a 
n 



r 



Q 
III 

a 
lii 

01 

a 

fi 

III 

CI 



> 
o 

o 




1000 



10 



T 1—1 

o 
o 



o 



1.1x10 
1.0x10-'° 
9.0x1 0-"" 



NtHse only (84S MHz) 

irK?in=»6,3E-ll 

standard deviadon = 3.7B-12 




200 



400 600 

Time index 



1000 



Fig- 10. 



